Change of quasiparticle dispersion in crossing Tc in the underdoped cuprates 
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One of the most remarkable properties of the high-temperature superconductors is a pseudogap 
regime appearing in the underdoped cuprates above the superconducting transition temperature Tc- 
The pseudogap continously develops out of the superconducting gap. In this paper, we demonstrate 
by means of a detailed comparison between theory and experiment that the characteristic change of 
quasiparticle dispersion in crossing Tc in the underdoped cuprates can be understood as being due 
to phase fluctuations of the superconducting order parameter. In particular, we show that within a 
phase fluctuation model the characteristic back-turning BCS bands disappear above Tc whereas the 
gap remains open. Furthermore, the pseudogap rather has a U-shape instead of the characteristic 
V-shape of a d^2_y2-wave pairing symmetry and starts closing from the nodal k = , ^) directions. 



whereas it rather fills in at the anti-nodal k 
fluctuation scenario. 



(tt, 0) regions, yielding further support to the phase 
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I. INTRODUCTION 

More than 15 years after the discovery of the high- 
Tc superconductors, the mechanisms leading to their un- 
usual properties are stiU under debate. Especially the 
pseudogap phase, which appears in various experiments 
below a characteristic temperature T* in the underdoped 
region of the phase diagram as a reduction of spectral 
weight^"^, might be a key to a better understanding of the 
high-Tc superconducting (SC) cuprates. In 1995 Emery 
and Kivelson^ proposed, that the proximity to the Mott 
insulating phase implies a strongly reduced phase stiff- 
ness J ^ compared to the usual BCS case. This 
causes the phase ordering temperature ^ J to be 
much lower than the mean-field pair-binding tempera- 
ture T^^^ . Taking this idea one step further'^ implies 
that at least part of the pseudogap behavior might be due 
to a kind of "pre-formed" Cooper pairs which form at a 
temperature T* = well above the actual SC transi- 
tion temperature Tc = T^, where phase coherence among 
these pairs finally sets in. This phase fluctuation scenario 
also explains quite natural the strongly enhanced Nernst 
signal above Tc in the underdoped cuprates-^-. 

In previous work, we have already shown that indeed 
a two-dimensional BCS-like Hamiltonian with a €1^2 _y2- 
wave gap and phase fluctuations, which were treated by 
a Monte-Carlo simulation of an XY model, yields re- 
sults which compare very well with scanning tunneling 
measurements over a wide temperature range^^. Fur- 
thermore, this phenomenological phase fluctuation model 
was also able to explain the possible "violation" of the in- 
plane optical integral in underdoped Bi2Sr2CaCu208+5 
(Bi2212)^ 

However, for the phase fluctuation description to be 
correct over a wide temperature range, one needs a mech- 



anism that produces "cheap" vortices, so that the only 
energy scale is the stiffness J and the dominating fluc- 
tuation channel is that of the phase of the SC order pa- 
rameter. Mechanisms that can lead to a small vortex 
core energy range from the more conventional picture of 
a granular superconductor, where the vortices arrange 
themselves to live in the insulating regions between the 
SC grains, up to the existence of a competing order that 
exists inside the vortex cores. As soon as the supercon- 
ducting order parameter is suppressed inside the vortex 
core, the system develops the competing order instead of 
going into a normal conducting paramagnetic state and 
thus has a much smaller vortex energy compared to a con- 
ventional BCS superconductor. Recently, it was shown&, 
that vortices with staggered-flux core can provide a way 
to understand the low vortex energy over a wide temper- 
ature range above Tc- In all cases, the small phase stiff- 
ness and the low vortex core energy have the same origin, 
which is the proximity to the Mott-insulating state. 

In this paper, we present theoretical results on the 
quasiparticle dispersion, which — when compared with 
experimental data — give a clear fingerprint towards 
a possible phase fluctuation scenario for the origin of 
the pseudogap. Earlier angle resolved photoemission 
(ARPES) results have shown deviations from the sim- 
ple BCS dx2-y2-'wave form of the SC gap in underdoped 
Bi22122i&, which might be due to a change in the pair- 
ing interaction in the proximity of the AF insulating 
phase. By analyzing the temperature dependence of the 
quasiparticle (QP) dispersion, we want to show, that the 
change of the QP dispersion in crossing Tc from the SC 
to the pseudogap region can be understood quite natu- 
rally by the assumption that the pseudogap is caused by 
phase fluctuations of the SC gap. Moreover, the phase 
fluctuation scenario also explains the deviations from the 



2 



simple rf2,2_j^2-wave form of SC and pseudogap^ in the 
underdoped cuprates. Using the ARPES with tunable 
excitation photon energy we disentangle bilayer splitting 
related effects and determine the true dispersion and the 
leading edge gap (LEG) function corresponding to the 
bonding band in the pseudogap regime of Pb-Bi2212. 

Since below Tc the QP excitations are perfectly BCS- 
like unless in the extremely underdoped region^", it is 
tempting to start from the BCS ground state and see how 
it is destroyed by including phase fluctuationsiitiSiiiii. 

II. MODEL AND CALCULATIONS 

In the following we use a phenomenological phase fluc- 
tuation model which has already been shown to success- 
fully account for the pseudogap observed in tunneling 
experiments'^ and which was also able to explain the pos- 
sible "violation" of the in-plane optical integral in under- 
doped Bi2213^. We consider the Hamiltonian 

H = Ho + Hi, (1) 

where Hq is the usual tight-binding Hamiltonian of non- 
interacting electrons on a two-dimensional (2D) square 
lattice 

^ = -^ E (4.^j.+4.^?J-^E-?.- (2) 

Here, cl (cj^.) creates (annihilates) an electron of spin a 

on the i site of the 2D square lattice and — cl c? 
is the number operator, t denotes an effective nearest- 
neighbor hopping-term and is the chemical potential. 
The angles (• • • ) indicate sums over nearest- neighbor 
sites of the 2D square lattice. 

The second part of the Hamiltonian }1\ contains a 
BCS-like d-wave interaction, which is given by 

i/, = -g^(A,^(At.) + At_(A,^)), (3) 

with b connecting nearest-neighbor sites. The coupling 
constant g stands for the strength of the effective next- 
neighbor (i2;2_y2-wave pairing-interaction. The origin of 
this pairing interaction is unimportant for the further 
calculation. It can be either of pure electronic origin, 
like spin fluctuations, or phonon mediated. The only 
important thing is, that there exists an effective pairing 
interaction, that produces a finite local d^2_y2-wave gap 
as one goes below a certain temperature T*. In contrast 
to conventional BCS theory, we consider the pairing-field 
amplitude not as a constant real number, but rather as 
a complex number 

(A1-) = ^(ct cl - -ct ct - ) = Ae**?^ (4) 

with a constant magnitude A and a fluctuating bond- 
phase field ^j^. In order to get a description, where the 



center of mass phases of the Cooper pairs are the only 
relevant degrees of freedomi^, we write the d^2_y2-wave 
bond-phase field in the following way 

{((fi^ + (fi^ , /2 for (5 in x-direction , , 

((pj -I- ip-'j^g)/2 + TT for 6 in y-direction, 

where is the center of mass phase of a Cooper pair 
localized at lattice site i. 

In order to account for the proximity to the Mott insu- 
lating state and thus the low superfluid density, we per- 
form a quenched average over all possible phase configu- 
rations with the statistical weight given by the classical 
XY free energy 

F[Vi] = - cos (if, -if j), (6) 

(ij) 

where the phase stiffness J determines the Berezinskii- 
Kosterlitz-Thouless transition temperature Tbkt to a 
quasi phase ordered state which we take as Tc- The XY 
free energy is defined on a coarse-grained lattice with the 
scale of the lattice spacing given by the pair coherence 
length^^ ^0 ^ . Now, the underdoped cuprates are 
in an intermediate coupling regime between large BCS 
mean-field pairs and tightly-bound BEC pairs^^'^^, with 
the pair-size coherence length given by 3 to 4 times 
the basic Cu-Cu lattice spacing. For a typical 36 x 36 
fcrmionic lattice, which is numerical feasible, we would 
only have a 9 x 9 phase lattice on top of it. This would not 
allow for any proper temperature scaling of the phase cor- 
relation length ^ (r) and obscure the Kosterlitz-Thouless 
transition. Therefore we have chosen to set A^c = 1.0 1. 
This yields i^o ^ 1 and allows the Monte Carlo phase 
simulation to be carried out on the same L x L lat- 
tice that is used for the diagonalization of the fermionic 
Hamiltonian^. In addition, the choice of Age — 1-0 i au- 
tomatically introduces the important short distance cut- 
off. Finally we set « \T*, where we had the STM 
experiments of Ref. in mind. 



III. EXPERIMENTAL DETAILS 

The ARPES experiments were carried out using angle- 
multiplexing electron energy analyzers. Spectra were 
recorded either with hi/ =21.218 eV photons from a He 
source or using radiation from the U125/1-PGM beam- 
line at the BESSY synchrotron radiation facility. The 
total energy resolution was set to 17 meV (FWHM) at 
hi/ = 38 eV. The angular resolution was kept below 0.2° 
both along and perpendicular to the analyzer entrance 
slit. Data shown in Fig. 0] were taken with 0.2° x 0.3° 
angular resolution. The data were collected on two sim- 
ilar underdoped, modulation-free single crystals of Pb- 
Bi2212 (Tc=77K). 
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FIG. 1: Spectral weight A{k,u}) in the pseudogap state (T = 
2.0 Tc, top) and in the superconducting state slightly below 
Tc {T = 0.75 Tc) calculated from the phase fluctuation model. 
For comparison we also show the spectral weight for the phase 
coherent BCS limit (bottom). 



IV. DISCUSSION OF RESULTS 
A. Dispersion 

Fig. n] shows the quasiparticle dispersion calculated 
from our phenomenological phase fluctuation model for 
10% doping {{n) = 0.9). The spectral weight is plot- 
ted along the (0, 0) (tt, 0) — > (tt, tt) direction through 
the Brillouin zone (BZ). The free dispersion would cross 
the Fermi surface close to the (tt, 0)-point. One can 
clearly see that the characteristic (back-turning) Bogoli- 
ubov quasiparticle band disappears in the pseudogap 
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FIG. 2: a) Superconducting state. Energy distribution of the 
photoemission intensity along the direction shown as a red 
arrow on the sketch below. The BCS-like dispersion is clearly 
observed for the bonding band, b) Pseudogap state. No more 
bending back of the dispersion is observed. Instead, spectral 
weight fades upon approaching the Fermi level. 



state above Tc- Instead, one obtains a sharp quasipar- 
ticle dispersion which runs straight towards the Fermi 
energy and than fades out at a distance of the order of 
the superconducting (SC) gap Asc- This is in complete 
agreement with the experimentally observed dispersion 
in underdoped Bi2212 which is shown in Fig. |2| 

The angle resolved data presented in Fig. |21 provide 
an insight into how the pseudogap is actually created 
near kf. In the superconducting state a characteristic 
BSC-like back-dispersion is easily seen. This clarity is 
achieved by the careful choice of the excitation photon 
energy. Exactly near hi^=38 cV the emission probabil- 
ity for the bonding band is much higher than for the 
antibonding ban d^^i^^ and bilayer-related complications 
are thus avoided. Above Tc in the pseudogap state the 
characteristic BCS behavior is replaced by the straight 
dispersion and strong depletion of the spectral weight 
towards Ef, which, as will be shown below, still leaves 
the energy gap in the spectrum. 

Furthermore, Fig. shows that the sharp quasiparti- 
cle features close to the (tt, 0)-point are getting lost above 
Tc within the phase fluctuation model. The sharp coher- 
ent (tt, 0)-peaks dissappear and broad incoherent weight 
EUs-in the gap. Exactly this behavior was observed be- 
fore in photoemission studies of the pseudogapi^iSfliSLSSiS^ 
and is also responsible for the characteristic temperature 
dependence of the scanning tunneling gap in the under- 
doped cuprateai^, where the pseudogap fills-in instead 
of closing. Interestingly, not only SC fluctuation^, but 
also staggered flux fluctuation^ can lead to this tem- 
perature dependence of the (tt, 0)-photoemission-peak. 
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The disappearance of the characteristic BCS bands 
above within the phase fluctuation picture can be un- 
derstood by the fact, that the BCS wave- function is a 
coherent superposition of wave functions with different 
number of electron pairs^^ 

\'^BCs)^X{{uk+Vkcl^cl^^)\M=Y,^N\^N), (7) 

k N 

where \^ n) is a N-particle wavefunction. The quantum 
mechanical uncertainty in the particle number is given 
by: 

{m^ = AY,ulvl. (8) 

k 

Now vf, — \ — u\ is the momentum distribution func- 
tion for T = and the weight of a quasiparticle peak at 
momentum k is given by vj, (w^) for E < Ej {E > Ef). 

In the normal metallic state with A = 0, one gets a 
sharp cut-off in (u^) at the Fermi- wavevector k = 
kf so that (AiV)^ = 0. In the BCS-superconducting 
state, however, {uD are finite also beyond the Fermi- 
wavevector kf, which means that also forfc > kf (k < kf) 
one gets spectral weight at E < Ef {E > Ef). This pro- 
duces the characteristic BCS band structure, with bands 
approaching Ef from below (above) and then turning 
back to higher binding (quasiparticle) energies. 

Now what happens if one introduces an arbitrary phase 
factor into the BCS wave function2& 

l*^>=n(l"'=l + l^'=|e'^4TcUi)l</'o). (9) 

k 

Integrating over all possible phases yieldsSS 
l*w>- r d^e-^^'^/'l[{\uk\ + \vk\e^''ci^cU^)\M 

Jo 

(10) 

This means, that one projects into an exact particle- 
number eigenstate by making the relative phase of the 
Cooper-pairs completely uncertain. Eq. ||SJ) is a special 
case of the general uncertainty relation between phase 
and particle number 

ANAip>l. (11) 

The above described behavior corresponds to what is 
happening in the phase-fluctuation model as a function 
of temperature. Starting from a phase coherent state at 
T = with Aif — 0, the particle number is completely 
uncertain with AN given by Eq.|Sl With increasing tem- 
perature, one gradually projects into a state with exact 
particle number N. In the temperature range, where the 
phases are completely uncorrelated (^ ^0)7 one then 
obtains AA'' = 0, and the back-turning BCS-bands must 
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FIG. 3: Gap function A(fc) in the pseudogap state (T = 
2.0 Tc, top) and in the superconducting state slightly below 
Tc {T — 0.75 Tc) calculated from the phase fluctuation model. 
For comparison we also show the BCS gap function (bottom) . 



completely dissappear (loose weight for k > kf, as seen 
in Figs. Hand 121). At finite temperatures, this situation 
corresponds to a classical grand canonical average over 
ensembles with different number of particles, where each 
state has a well defined particle number and is no longer 
a coherent quantum-mechanical superposition of states 
with different number of particles. Thus, we obtain a 
crossover from a BCS-like phase-ordered band structure 
to a completely new phase-disordered pseudo-gaped band 
structure. 
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FIG. 4: Values of the leading edge (pseudo-)gap (LEG) as 
a function of the Fermi surface angle within the quadrant of 
the Brillouin zone (see lower panel) given with respect to the 
binding energy of the leading edge of a nodal energy distribu- 
tion curve. The curve joining the low-gap extremity of these 
data points would represent the fc-dependence of the pseudo- 
gap. For details see Ref. |^ 

B. Superconducting gap and pseudogap 

Next we want to elucidate the effect of phase fluctu- 
ations on the fc-dependence of the quasiparticle pairing 
gap. Therefore, we have plotted in Fig.Othe quasiparti- 
cle dispersion obtained from MC simulations of the phase 
fluctuation model along the Fermi surface of the free dis- 
persion e{k) at half-filling ((n) — 1.0). This gives us ef- 
fectively the gap function A (A;). As can be seen in Fig.O 
below Tc one obtains the characteristic F-shape of a gap 
with dj;2_j^2 pairing symmetry. As the temperature is 
raised, the quasiparticles peaks are getting broader. In 
the pseudogap state above Tc, the spectral weight is get- 
ting rather incoherent close the the k — (tt, 0) as was 
pointed out before. However, close to the nodal point 
of the gap function one still obtains a sharp quasiparti- 
cle dispersion. There, one can clearly see spectral weight 
shifting to lower binding energies which produces an ex- 
tended gapless region in the pseudogap state close to fc = 
instead of the nodal point in the superconducting 
state below Tc- This behavior is in complete agreement 
with photoemission experiments^'^"^ which show that the 
pseudogap starts closing from k — (f , f ) where one ob- 



tains a finite Fermi-arc but rather Ells in at fc = (tt, 0) 
exactly as in Fig.|31(top). 

Furthermore, the pseudogap A(fc) obtained from phase 
fluctuations of the local dx^-y^ pairing-amplitude rather 
has a U-like shape (see Fig. O top) than the charac- 
teristic V-shape of a BCS dx2_y2-geij>. For comparison, 
Fig. 0] shows the experimentally observed pseudogap in 
underdoped Bi2212. One can clearly see that the exper- 
imentally determined pseudogap has exactly the U-like 
form that we have obtained from the phase fluctuation 
model. This deviation from the pure da;2_j,2-wave form 
was also observed in the superconducting state of very 
underdoped cuprates and interpreted as higher harmonic 
contributions to the pairing function^ '*. However, these 
experimental results just might indicate the possible rel- 
evance of quantum phase fluctuations in this region of 
the phase diagram. Our results on the effects of phase 
fluctuations on the form of the pairing gap could also be 
of some relevance for electron doped cupratcs^^ where a 
possible crossover from a dx2_y2 (or anisotropic s-wave) 
to a pure s-wave symmetry of the superconducting gap 
as a function of electron doping was observed^SiSS. 

V. SUMMARY AND CONCLUSION 

In conclusion, we have elaborated the important role 
that phase fluctuation effects might play in the under- 
doped cuprates. With a detailed comparison between 
theory and experiment we were able to show how phase 
fluctuations influence the quasiparticle spectra. In par- 
ticular the disappearance of the BCS-Bogoliubov quasi- 
particle band at Tc and the change from a more V-like 
superconducting gap to a rather U-like pseudogap above 
Tc can be explained in a consistent way by assuming that 
the low energy pseudogap in the underdoped cuprates is 
due to phase fluctuations of a local dj,2_j,2-wave pairing 
gap with fixed magnitude. Furthermore, phase fluctua- 
tions can explain why the pseudogap starts closing from 
the nodal points, whereas it rather fills in along the anti- 
nodal directions. 
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